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DISCRETE  FOURIER  TRANSFORMS 
FOR  SYMMETRIC  OR  ANTISYMMETRIC  REAL  SEQUENCES 
WITH  THE  POINT  OF  SYMMETRY  OR  ANTISYMMETRY 
MIDWAY  BETWEEN  TWO  CELLS 


1.  Introduction 

In  this  paper  we  present  a method  for  computing  the  discrete  Fourier  transform  (DFT)  of 
real  sequences  (>(>))  of  length  AN  which  satisfy  the  symmetric  or  antisymmetric  boundary 
conditions 


y(4\  -j  - l)  -sy(J),  O^j^AN  -1.  s - ±1. 


(1.1) 


Such  sequences  appear  in  several  important  classes  of  problems  in  which  certain  considerations 
(physical  boundaries,  reflection  points,  etc.)  dictate  that  the  symmetry  or  antisymmetry  points 
be  at  the  midpoint  between  cells.  For  example,  such  sequences  are  involved  in  solving 
Poisson's  equation  V 2<f>  - p for  the  gravitational  potential  <t>K  or  the  electrostatic  potential  4>,,, 
given  the  mass  or  charge  density  when  the  source  terms  are  specified  half  a cell  away  from  the 
points  of  symmetry  or  antisymmetry.  The  conditions  (1.1)  reflect  this  property,  supposing 
symmetry  or  antisymmetry  of  the  sequence  {>(>))  at  points  halfway  between  >(-1)  and 
>(0),  y(2N  -1)  and  y(2N),  and  y(4N  —1)  and  y(AN).  This  contrasts  with  situations  for 
which  the  usual  sine-cosine  transforms  are  applicable,  in  which  >(0),  y(2N),  and  y(4N)  are 
the  points  of  symmetry  or  antisymmetry  and  have  the  symmetric  or  antisymmetric  boundary 
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y(4N  — j)  - sy(j).  O^j^AN.s  — ±1. 


(1.2) 


The  approach  presented  here  is  more  appropriate  for  efficient  calculation  of  DFTs  of  real  se- 
quences satisfying  (1.1)  than  the  sine-cosine  transforms  of  [1]. 

Our  technique  for  computing  the  DFT 


4A-1 


Y(k)  - £ yij)v(jk),  0<fc<4A(-l. 

)- o 


(1.3) 


where  w(q)  - exp  ( -2iriq/4N),  approximately  halves  the  number  of  operations  which  are  re- 
quired by  algorithms  in  which  the  added  conditions  of  symmetry  or  antisymmetry  are  not 
taken  into  account. 

2.  Derivation  of  the  Algorithm 

The  DFT  of  any  real  sequence  \y(j))  consisting  of  an  even  number  of  points 
2 M ( M - 27V  in  our  case)  may  be  calculated  by  performing  a DFT  on  M complex  points,  plus 
a little  overhead.  Specifically,  if 

x(j)  ~y(2j)  + iy(2J  + 1), 


(2.1) 


and 


M-l 


X(k)  - £ x(j)  exp(  -2niJk/M).  0<tc<Af-l. 

y-o 

then  it  is  a straightforward  exercise  to  verify  that 

F<*)  - 11  -iw(k)]X(k)/2  + [1  +iw(k))X*(M  — k)/2, 

Y(2M  - k)  - Y*(k),  0<  A/. 


(2.2) 


(2.3) 


where  X(M)  - 3f(0)  and  p*  denotes  the  complex  conjugate  of  p.  For  inversion,  the 
coefficients  Mf (Ac))  can  be  constructed  from  the  coefficients  ( Y(k))  and  are  given  by 
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X(k)  - (1  + /w(  -k)]Y(k)/2  + 11  -iw(-k))Y'(M  -k)/2,  0^j<M-l.  (2.4) 

The  boundary  values  on  Che  sequence  {*(./))  induced  by  (1.1)  are  given  by 

xW-y-1)  -/»u  o<y<  A/  -1.  (2.5) 

Therefore,  for  0 < k < M - 1 , 

**(*)  - “£  x'(j)w(-2jk)  - -is  *£  x(A/ -1)h-(-2j<c) 
y-o  y -o 

is  *£  jr(yM  -2k(M  -J -1)1 isw(2k)X(k).  (2.6) 

j- o 

From  (2.5)  it  follows  that  (2.2)  can  be  written  as  a sum  of  two  DFTs,  each  of  length  N: 

N- 1 N-l 

X(k)  - £ x(2j)w(4Jk)  + w(2k)  £ x(2y  + 1 )w(4jk) 

y •*>  )*> 

- "£  x(2j)w(4jk)  + isw(-2k)  *£  x'(2J)w(-4jk). 

> y-o  y -o 

Therefore,  letting 

B(*)  - "£  x(27)w(4y*),  (2.7) 

y-o 

and  noting  that  + *)  - 5(A),  we  have 
*(*)  - fl(*)  + /sw(  -2k)B*(k). 

X(N  + k ) -fl(*)  ~/sw(  -2k)B*(k). 

X(2N—k)  — B(N  — k)  + isw(2k)B'(N  — k).  OKk^N-l.  (2.8) 

• Hence,  all  4N  Fourier  coefficients  | Y(k))  can  be  found  from  a transform  of  the  N complex 

points  U(2y):  0 <J < N —\)  -ly(4»  + iy(4J  + 1 ):  0 < N - 1)  . 

\ 

Using  the  expressions  in  (2J)  to  evaluate  the  quantities  ( K(*)l  in  (2.3),  we  obtain 
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Y(k)  - w(  -k/2)  [Q(k)  + sQ'(k)]/2,  0<k<M-l.  (2.9) 

where 

Q(k)  - [w(k/2)  -iw(3k/2)]B(k)  + \w(k/2)  + iwttk/2 ) )B'(N -k). 

Q(N  + k)  - M<W  + *)/2}  — /wl3  (N  + k)/2])B(k) 

+ M(A'  + *)/2)  +iw[3(N  + k)/2))B'(N -k).  O^k^N-l.  (2.10) 

We  should  note  that  the  quantities  Q(k)  +sQ*(k)  in  (2.9)  are  either  real  (s  — 1)  or 
purely  imaginary  (s  — -1).  Because  of  the  symmetry  conditions  (1.1)  there  are  only  2 N linear- 
ly independent  coefficients.  The  N complex  coefficients  ffl(A))  are  one  such  set,  while  the  2 N 
real  or  purely  imaginary  coefficients  | Y(k))  defined  by 

Y(k)  - l(?(*)  +sQ*(*)]/2,  O^k^M-l.  (2. II) 

are  another.  Notice  that  the  coefficients  { K(A)|  and  |f(A)|  differ  only  by  the  phase  factor 
w(  -k/2).  It  is  a straightforward  matter  to  show  that  (K(A)l  is  the  set  of  Fourier  coefficients 
for  the  pure  sine  or  cosine  expansion  of  the  original  sequence,  if  a coordinate  system  offset  by 
half  a cell  is  used.  This  does  not  allow  the  coefficients  to  be  readily  calculated  by  the  methods 
of  111,  since  in  this  shifted  representation  the  original  sequence  points  are  now  specified  one- 
half  cell  from  the  expansion  points.  For  inversion,  the  coefficients  1 B(k))  can  be  expressed  in 
terms  of  the  quantities  ( K(A))  as  follows: 

4 B(k)  - lw(-k/2)  + iw(  -2k/2)]Y(k) 

+ MOV  - k)/2 1 + iw  13  OV  - k)/2 1|  Y*(N  - k) 

+ (w{-(W  + A)/2l  + M ~3(iV  + k)/2))Y(N  + k) 

- liw(-k/2)  + w(-3k/2))Y'(2N  - k).  0^k^N-\.  (212) 

Of  course,  the  decision  to  utilize  a particular  one  of  the  sets  |/f(/c)),  I K(A)),  or  1 K(A)}  depends 
upon  the  particular  problem  being  treated.  If  computing  speed  is  the  primary  consideration, 
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then  working  directly  with  the  quantities  {£(*))  is  the  best  approach.  However,  for  involved 
manipulations  of  the  transform  coefficients,  {K(*»  and  (TU)}  are  more  easily  interpreted 
physically,  and  use  of  one  of  these  sets  may  be  more  appropriate. 


Frequently  the  DFT  is  performed  to  facilitate  "smoothing"  the  data,  or  "filtering"  it  in 
some  way.  This  is  often  accomplished  by  multiplying  the  Fourier  coefficients  {K(/c)|  by  "form 
factors"  FUc),  0 ^ k^4N  -1,  which  are  all  real  or  all  purely  imaginary,  and  then  performing 
the  inverse  DFT 

, 4N-1 

Z 0)  -yr:  £ Hk)nk)w(-jk).  0<X4Af-l. 

*-4) 

In  the  present  problem,  we  impose  the  restrictions 

F{k)  - F'(*N  - k).  1<*<4AT-1.  (2.13) 


We  write 


F*(k)  ~tF(k),  0<*<4JV-1. 

where  t - +1  or  -1  according  as  each  form  factor  is  real  or  purely  imaginary,  and  F(4JV)  - 
F( 0).  Using  (2.3)  it  is  easily  verified  that 


z*(j)  -z(J).  0<J<AN  -1,  (2.14) 

and  hence  the  "smoothed"  coefficients  (z(y))  are  real.  Using  (2.3)  and  (2.13)  it  may  be  directly 
shown  that 

z(AN  — J — \)  — stz(J),  0^j<AN-\.  (2.15) 


Thus,  the  sequence  (z(y)l  is  a member  of  the  symmetry  class  defined  by  (1.1),  so  from 
(2.5)  and  (2.7)  we  may  conclude  the  existence  of  a sequence  |2>'(*)}  such  that 


*'(*>  -4;  l'  u(lj)w{AJk).  0*k*N-l. 
N )• 0 
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where 

u(M  -j  - 1 ) - istu*(j), 

u(j)  — z(2j)  + iz(2j  + 1 ),  0 < j < A/ — 1. 

Consequently,  the  "smoothed"  data  points  \z(J))  can  be  determined  from  an  inverse  transform 
of  length  N applied  to  the  coefficients  (B'(/c)).  Specifically, 

u(2j)  B'(k)  w(  —4jk), 

N k *0 

It  remains  to  express  the  quantities  {/?'(*))  in  terms  of  the  coefficients  {B(/c)|.  From 
(2,12)  we  obtain 

AB(k)  - [w(-k/ 2)  + iw(-3k/2)]Y'(k) 

+ { w(  (A(  - k)/2 1 + iw[UN  - k)/2]\Y’\ 

+ {w[-(N  + k)/2]  + /w(  —3  (N  + k)/2  ])  Y'(N 

- [/>(  -k/2)  + w(  -3k/2)]Y'*(2N  -k),  0<ft<Af-I,  (2.16) 

where 

Y'(k)  - F(k)Y(k),  0<£*<M-1.  (2.17) 

Substituting  (2.17)  into  (2.16),  it  is  a bit  tedious  but  not  difficult  to  verify  that 

8 B (k)  - akB(k)  + pkB*(k)  + ykB(N  - k)  + 6kB*(N  - k). 

where 

a*  -2  [Fik)  +F*(2N-k)  + F(N  + k)  + F*(/V-A)] 

-ilF(k)  -F'(2N -k))lw(k)  — w(  —k)  ] 

— IFIN  + k)  -F*{N  -k)]lw(k)  + *(-*)]; 
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(ik  - /sw(  -2k)[2\F(lc)  + F*(2N -k)  — F(N  + k)  -F*(N-k)  ] 

-HF(k)  — F*(2N  -*)]  Mfc)  w ( -fc)] 

+ [F(N  + k)  — F*{N -k))lw(k)  + w(  -fc)]); 

yk  - isw(  —2k){  -HF(k)  — F*(2N  — At)Hh»(*)  + w(  -*)) 

+ [F(N  + k)  -F'(N  -k))[w(k)  -tv(  -*))}; 

hk  - i[F(k ) -F'(2N  -k)][w(k)  +w(-k)) 

+ [m  + *>  -F*(N -k)][w(k)  -w(  -k)\. 

Consequently,  calculation  of  the  quantities  (/(*)}  or  { Y(k))  is  not  necessary.  All  that  is 
required  is  to  compute  the  transform  {£(&)},  obtain  the  quantities  l0'(k)|  using  the  stored 
constants  |a  J,  {0*},  {>-*),  («*},  and  compute  the  inverse  transform  {u(2J)\. 

3.  Conclusions 

The  4Wreal  data  points  [y(j)\  0 AN  -1)  satisfying  the  symmetric  or  antisymmetric 
boundary  conditions  (1.1)  may  be  "processed"  (transformed,  smoothed,  and  inverse 
transformed)  using  two  DFTs  of  length  N.  This  is  a savings  of  approximately  a factor  of  two  in 
the  number  of  arithmetic  operations  over  the  computation  required  by  DFT  algorithms  in 
which  symmetry  or  antisymmetry  is  not  taken  into  account.  In  two  and  three  dimensional  ap- 
plications, we  obtain  savings  factors  of  four  and  eight,  respectively. 
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